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Abstract 

A simple experimental protocol applying a quantitative ultrasound (QUS) pulse-echo technique was used to measure the 
acoustic parameters of healthy femoral diaphyses of Wistar rats in vivo. Five quantitative parameters [apparent integrated 
backscatter (AIB), frequency slope of apparent backscatter (FSAB), time slope of apparent backscatter (TSAB), integrated 
reflection coefficient (IRC), and frequency slope of integrated reflection (FSIR)] were calculated using the echoes from cortical 
and trabecular bone in the femurs of 14 Wistar rats. Signal acquisition was performed three times in each rat, with the 
ultrasound signal acquired along the femur's central region from three positions 1 mm apart from each other. The parameters 
estimated for the three positions were averaged to represent the femur diaphysis. The results showed that AIB, FSAB, TSAB, 
and IRC values were statistically similar, but the FSIR values from Experiments 1 and 3 were different. Furthermore, Pearson's 
correlation coefficient showed, in general, strong correlations among the parameters. The proposed protocol and calculated 
parameters demonstrated the potential to characterize the femur diaphysis of rats in vivo. The results are relevant because 
rats have a bone structure very similar to humans, and thus are an important step toward preclinical trials and subsequent 
application of QUS in humans. 
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Introduction 

Bone is formed by specialized connective tissue; its 
extracellular matrix is calcified but maintains a degree of 
elasticity (1,2). Despite its hardness and resilience owing 
to an association between collagen and hydroxyapatite 
crystals (2,3), bone injury is still a recurrent health 
condition. The process of fracture healing is complex 
(3,4), involving cellular proliferation and differentiation, 
chemotaxis, and synthesis of extracellular matrix (3), with 
several stages of repair with a well-defined temporal and 
spatial sequence (5,6). Severe complications frequently 
occur during the process of bone repair (7-9), such as 
delayed union, malunion, or nonunion (pseudarthrosis) 
(10), resulting in negative consequences for patients and 
increased costs to healthcare systems (7,8), and justifying 
efforts to develop new diagnostic tools to track the bone 
healing process. 

Protopappas et al. (11), using computational simulations 
and experiments (4), found that material properties and 
geometrical features change, influencing the propagation of 



ultrasound waves along the cortex of long bones. Dodd et al. 
(12) also pointed out that fracture gaps can promote loss 
of ultrasound energy. Hakulinen et al. (13) studied the 
measurement of bone density and mechanical properties by 
ultrasonic propagation parameters, reporting results that 
indicated a potential for diagnosing osteoporosis. 

Application of quantitative ultrasound (QUS) to mea- 
sure acoustic parameters of bone structure (2,5,6,8,14, 
15) can enable the detection of changes and information 
about deviations from the normal condition caused by 
diseases (14,16-18) such as osteoporosis (13,19,20), 
osteomyelitis (10), and osteoarthritis. QUS can potentially 
minimize subjectivity in diagnoses made by conventional 
imaging methods. It has the advantage of low operat- 
ing costs (13), ease in equipment handling (16,21), and 
employing nonionizing radiation (22), unlike conventio- 
nal X-ray and computed tomography (23). Although 
the literature reports innumerable uses of QUS for the 
characterization of soft biological tissue (14,24), clinical 
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application is still modest (21) because of the lack of 
reproducibility of the method. Nevertheless, QUS is 
emerging as a promising modality for bone character- 
ization (21). Methods to characterize bone have been 
developed using animal models, such as cattle (19,25) 
and sheep (1 1 ,26), and have been subsequently explored 
in humans (15,18,20). Among animal bone models, few 
studies (27) have explored the use of QUS to quantify 
reflection and backscattering from the bones of Wistar 
rat models in vivo, which are the most similar to human 
bones. 

The present study aimed to develop a QUS method 
using a pulse-echo technique to extract five parameters: 
apparent integrated backscatter (AIB), frequency slope 
of apparent backscatter (FSAB), time slope of apparent 
backscatter (TSAB), integrated reflection coefficient 
(IRC), and frequency slope of integrated reflection (FSIR) 
by processing ultrasonic backscattered and reflection 
echoes to characterize cortical and trabecular bone from 
healthy bone diaphyses of Wistar rats in vivo. Although 
there is no adequate standard characterization method for 
diaphyseal bone in the literature, there exist two previous 
conference papers published by our group on this subject 
(27,28). To our knowledge, no other similar in vivo studies 
in rats have been published. 

Material and Methods 

The research was approved by the Ethics Committee 
for the Use of Laboratory Animals in Research of the 
Faculdade de Medicina, Universidade Federal do Rio de 
Janeiro. The animals were housed in accordance with 
the Guidelines for Care and Use of Animals in Research. 
The sample consisted of seven 3-month-old Wistar rats 
(Rattus norvegicus albinus) weighing 225 ± 25 g, previously 
anesthetized, and with their two hind legs shaved. 



by maximizing the echo signals (displayed on an 
oscilloscope screen). The central point position corre- 
sponded to the maximum echo and two other points, 
displaced 1 mm to each side, were also chosen. Signals 
from the three points were recorded. 

A reference signal from a polished steel plate (1 cm 
thick) was acquired by pointing the waveguide to the 
plate's front surface, which was covered with gel, but 
without the animal sample. 

Soft tissue thickness was measured with a B-mode 
scanning VEVO" 770 (VisualSonics, Inc., Canada) 
transducer at 30 MHz. The time position of the echo from 
the muscle/bone interface was identified in the radio- 
frequency (RF) signal, according to the transducer-bone 
distance calculated from the B-mode images. The 
experiment was carried out three times, with an interval 
of 10 days, and the experimental setup was rebuilt each 
time before measurements took place. 

Measurements of ultrasound parameters 

The transducer (model V326, Olympus® NDT, Inc., 
USA) with a nominal frequency of 5 MHz, 9.5 mm in 
diameter, and 69.3 mm in focal length, was driven by a 
pulse generator (model SR9000, Matec®, Inc., USA). The 
echoes were displayed and measured on an oscilloscope 
(model TDS 2024B, Tektronix 1 *, Inc., USA). 

The RF echo segment (Figure 1) used for character- 
ization of the bone diaphysis consisted of the following: i) 
the echo from the bone surface and ii) the backscatter- 
ing signal of the inner bone. The window containing this 
segment was determined by identifying the echo from the 
bone surface, which corresponded to the reference echo 
from the steel plate, and selecting a first rectangular 
window encompassing its limits using a 10% peak 
amplitude threshold (Figure 2). The segment of 4-u.s 
duration containing the backscattered signal from the 



Experimental setup 

The acquisition of ultrasound signals was performed 
by the same researcher three times in each animal at 
an interval of 10 days under the same environmental 
conditions and temperature (22.5±1.1°C). The animals 
were killed following the last measurement. 

The protocol for signal acquisition was as follows: a) 
Each animal was laid down in a lateral position with its 
hind limb relaxed and resting on a polished steel plate 
placed perpendicular to the ultrasonic beam axis, b) The 
transducer was held vertically by a stereotactic holder of 
a 2-\xvn resolution, with the beam focused on the femur 
diaphysis (region of interest), c) The coupling of transduc- 
er to the limb covered with a water-soluble gel was made 
with a waveguide, a glass tube filled with degassed water 
sealed by 10.5-um thick polyvinyl chloride, d) Adjustment 
of the transducer position was guided by palpation (the 
greater trochanter and lateral condyle of the femur were 
found to define the lateral middle third of the femur) and 
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Figure 1. Radiofrequency signal displaying the echo interfaces. 
A, Degassed water/PVC membrane interface; B, skin/fat/muscle 
interface; C, window from surface bone (reflection, continuous 
line box); D, window from inside bone (backscatter, dotted line 
box); £, surface bone/skin/fat/muscle interface (limb/gel); F, skin/ 
reflective plate interface. 
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Figure 2. A, Definition of the rectangular window over the 
reflected echo from the reference steel plate (extreme limits on 
1 0% of peak amplitude). B, Reference rectangular window placed 
over the echo from interface muscle/bone. The remaining echoes 
on the right with a 4-us duration are from backscattering inside 
the bone. 

inside of the bone (which has no interference from the 
muscle/bone interface echo) began from the end of this 
window. This time duration was chosen to ensure that the 
backscattering signal came from within the bone, which, in 
adult Wistar rats, has an estimated average diameter of 
3.12±0.1 mm. 

Five parameters were used to characterize the echoes 
from specimen and reference signals: AIB, FSAB, TSAB, 
IRC, and FSIR, and an algorithm was developed in 
Matlab® code (MathWorks, Inc., USA) to calculate these 
parameters. The apparent backscatter transfer function 
(ABTF; Equation 1) was obtained from the literature (20) 
where Preference and P sp ecimen are the power spectra of 
reference and specimen signals (rectangular window from 
bone surface), respectively. 

ABTF= 1 0 log 10 P spec , men (0 - 1 0 log 10 P re ference{f) (1) 

AIB, determined by integrating the curve of ABTF 
(Equation 2), expresses the average value of apparent 
backscattering in the frequency range of interest, while 
FSAB, the slope of the linear regression line obtained 
from the curve ABTF vs frequency, is the fraction of 
apparent backscattering corresponding to each fre- 
quency. Furthermore, TSAB, the slope of the linear 
regression of values of AIB as a function of time, is the 
variation of apparent backscattering as the wave propa- 
gates through the tissue. 

flugh 

[ABTF].df 

AIB= — (2) 

'high 'low 

Reflection transfer function (RTF; Equation 3) was 
calculated from the rectangular window with the reflection 
signal, which has a definition similar to ABTF, while IRC 
is the average value of reflection within a frequency 
bandwidth, which is obtained by integrating RTF. FSIR is 
the slope of the linear regression of IRC and is the 
apparent fraction of the reflection corresponding to each 
frequency, and calculated similar to FSAB (20). 




Figure 3. Layers of biological tissues (skin, fat, muscle, cortical 
bone, inside bone) and reflective plate. 

RTF= 1 0 log 10 P S pecimen(f) ~ 1 0 log 10 P re ference(f) (3) 

It must be noted that, before the estimation of 
parameters, all RF echoes must have their amplitudes 
corrected for loss in amplitude because of attenuation in 
the propagation path. By assuming that the medium is 
composed of four layers (skin, fat, muscle, and bone; 
Figure 3), the compensation factor A f can be written as 
follows (Equation 4): where s, f and m indices correspond 
to skin, fat, and muscle layers, respectively, while as, af 
and am denote their respective attenuation coefficients, 
and xs, xf, and xm indicate their respective thicknesses. 
Similarly, T s/f and 7>/ m are the transmission coefficients 
for the skin/fat and fat/muscle interfaces, respectively. 
The compensation factor takes into account the propaga- 
tion to and from the bone surface. 

a «2as.xST-2 Jlaf.xfT?- a 2am.xm /a\ 

H f— e 's/f e 'f/m e W 

Five parameters were calculated from the three signals 
for each femur, and the average of these values was 
considered as representative of each hind limb. Thus, 
each femur was characterized by five parameters. 

Quantitative computed tomography (QCT) 
acquisition 

QCT (75 kV, 145 mAs) was performed with Triumph II 
PET/SPECT/micro-CT equipment (Gamma Medica, Inc., 
Canada). Eight femurs were put in the scanning plate in 
a lateral position (the same one used for ultrasound 
acquisition). The tomographic images were processed 
with the Osirix software for bone density analysis (in 
Hounsfield units) of the femur diaphysis. These data were 
used as the gold standard for bone density. 

Statistical analysis 

The Kolmogorov-Smirnov and equal variance tests 
were used to test normality; when these tests failed, the 
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Table 1. Average values and standard deviation of each parameter in three experiments. 



Experiments 


AIB (dB) 


FSAB (dB/MHz) 


TSAB (dB/u.s) 


IRC (dB) 


FSIR (dB/MHz) 


1 


-44.46 ± 2.23 


-1.43 ± 0.28 


-1.31 ± 0.23 


-28.07 ±1.17 


-1.62 ± 0.08 


2 


-43.96 ± 1.09 


-1.39 ± 0.27 


-1.30 ± 0.16 


-28.13 ± 0.73 


-1.63 ± 0.05 


3 


-43.76 ± 1.93 


-1.31 ± 0.26 


-1.22 ± 0.18 


-28.46 ± 0.85 


-1.66 ± 0.06 



AIB: apparent integrated backscatter; FSAB: frequency slope of apparent backscatter; TSAB: time slope of apparent backscatter; IRC: 
integrated reflection coefficient; FSIR: frequency slope of integrated reflection. 



nonparametric Friedman test was used. To test the null 
hypothesis that all parameters belonged to the same 
population, and to assess the reproducibility of the 
parameters and method used, a statistical analysis was 
performed by one-way repeated measures analysis of 
variance (ANOVA) at a 5% level of significance. The tests 
were carried out using SigmaStat 3.5 (Systat Software, 
Inc., USA). Pearson's correlation coefficient was used to 
quantify the correlation between the reflection parameters 
(IRC and FSIR) and surface bone density (as measured 
by QCT). Pearson's correlation was also used to evaluate 
the level of association between the ultrasonic para- 
meters. 

Results 

The average values and standard deviations of 
parameters AIB, FSAB, TSAB, IRC, and FSIR for each 
experiment are shown in Table 1. 

Pearson's correlation coefficient showed a positive 
correlation between surface bone density and IRC (Figure 
4A) and FSIR (Figure 4B) for the three experiments. 

To determine the interaction between the parameters 
from backscatter (AIB, FSAB, and TSAB; Figure 5) and 
reflection (IRC and FSIR) of bones, Pearson's correlation 
was applied to parameters of Experiment 1 (Table 2), Ex- 
periment 2 (Table 3), and Experiment 3 (Table 4). The one- 
way repeated measures ANOVA was used independently 



to test the AIB, FSAB, TSAB, and FSIR data from the three 
experiments, which showed that AIB [F(2,41) = 1.02; 
P = 0.37], FSAB [F(2,41) = 1.92; P = 0.17], and TSAB 
[F(2,41) = 1.61; P = 0.22] belong to the same population. 
With regard to the FSIR parameter, the differences in mean 
values among the experiments were greater than those 
expected, thus there was a statistically significant difference 
[F(2,41) = 4.17; P = 0.03]. To isolate the experiment that 
differed from the others, a pairwise multiple comparison 
procedure (Holm-Sidak method) at a 5% level of signifi- 
cance was employed, which showed that the FSIR 
parameter of Experiment 1 was different from that of 
Experiment 3 (P = 0.01). 

The IRC parameter did not follow a normal distribu- 
tion, and thus we used the Friedman test at a 5% level of 
significance, which showed that this parameter belonged 
to the same population [F(2) = 4; P = 0.14]. Pearson's 
correlation coefficient showed, in general, strong correla- 
tions among the parameters. 

Discussion 

The literature on new protocols and tools for the 
characterization of bone by QUS is extensive (19- 
21,26,28), but, to date, there is no standardization of 
methods. In previous studies, several analyses have 
been conducted, including mathematical simulations 
(11), in vitro studies (12,13,25), and even animal models 
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Figure 4. Pearson's correlation between surface bone density, the integrated reflection coefficient (IRC; A) and frequency slope of 
integrated reflection (FSIR; S) for three experiments. The experiment showed positive correlation between surface bone density and 
the parameter IRC (r = 0.79, P = 0.019; r = 0.87, P = 0.005; r = 0.724, P = 0.042), in Experiments 1, 2, and 3, respectively, and the 
parameter FSIR (r = 0.78; P = 0.023; r=0.77; P = 0.026; R = 0.81; P = 0.016), in Experiments 1, 2, and 3, respectively. 
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Figure 5. Correlation between the backscatter parameters (AIB, FSAB, and TSAB) in Experiments 1-3 of the bone. AIB: apparent 
integrated backscatter; FSAB: frequency slope of apparent backscatter; TSAB: time slope of apparent backscatter. 



(26,28), with different degrees of success. Even the 
promising studies are noted to have difficulties in 
extrapolating their results for use in preclinical trials. 
The following are some of the reasons for this limitation. 
The bone specimens examined had important differ- 
ences, when compared with human bone (17,19,25). The 
simulation models were very simple, or the method 
cannot be applied in vivo; and most of all, there is an 
important lack of reproducibility. When dealing with 
human in vivo experiments, studies are carried out mostly 
in a transmission mode applied to extremities (calcaneus, 
phalange, and forearm) (21), and the same difficulty is 
encountered with regard to definition of protocols as well 
as reproducibility. 

QUS can help in minimizing the subjectivity of bone 
imaging and characterization and is more accurate for 
diagnosis and monitoring either the evolution of treatment 
or the course of metabolic diseases. In the present study, 
an in vivo experiment was proposed along with standardi- 
zation and consolidation of protocols as a step prior to 
clinical use in humans. The first proposal was the use of live 
adult Wistar rats that, in addition to their easy handling and 
resistance in research (9,29), have bone tissue character- 
istics more similar to human bones (29,30) than other 
animals, with the exception of primates. Furthermore, as 
rats play an important role in the evaluation of metabolic 
bone diseases (29) and pathophysiological conditions 
(29,30), the results are more likely to be similar to those 



expected in humans. Because rat femurs are small, it was 
necessary to use a 5-MHz transducer to provide satisfac- 
tory resolution (20). Additionally, the precise positioning of 
the ultrasonic beam at the three acquisition sites on the 
bone, chosen close enough to ensure minimal anatomical 
variations, was ensured by a high-precision stereotactic 
holder. The acquisition of in vivo signals is highly influenced 
by soft tissue attenuation and reflection; therefore, a com- 
pensation factor (A f ) to reduce measurement errors (31) 
was adopted. 

QCT provides accurate bone density measurements 
of the cortical bone that is used as a gold standard. In the 
three experiments, a positive correlation between the IRC 
and FSIR parameters and surface bone density was 
identified. According to the value correlation coefficient, 
we consider our method as adequate to characterize 
cortical bone density. 

The five parameters estimated were shown to be 
appropriate for characterizing the rat femurs in vivo. It is 
interesting to note that the first three parameters were 
described by Hoffmeister et al. (20) for in vitro character- 
ization of human cancellous bones. On the other hand, IRC 
was used for in vitro characterization of trabecular bovine 
bone (13,32) and in vitro trabecular human bone (31,33), 
and FSIR was developed by our group for characterization 
of rat bone. The AIB parameter can be used to characte- 
rize soft tissues (24). Hoffmeister et al. (20) used AIB 
(-40.9 ± 2.0 dB), FSAB (-1.7 ±0.5 dB/MHz), and TSAB 



Table 2. Pearson's correlation coefficient (r) of parameters and P for Experiment 1 . 





AIB (dB) 


FSAB (dB/MHz) 


TSAB (dB/u.s) 


IRC (dB) 


FSIR (dB/MHz) 


AIB (dB) 
FSAB (dB/MHz) 
TSAB (dB/us) 
IRC (dB) 


r = 0.92 c 
r = 0.98 d 
r = -0.87 b 


r = 0.92 c 

r = 0.89 b 
r = -0.93 c 


r = 0.98 d 
r = 0.89 b 

r = -0.86 b 


r = -0.87 b 
r = -0.93 c 
r = -0.86 b 


r = -0.83 a 
r = -0.86 b 
r = -0.79 a 
r = 0.85 a 



a P<0.001; b P<0.0001; C P<0.00001; d P<0.00000001 . AIB: apparent integrated backscatter; FSAB: frequency slope of apparent 
backscatter; TSAB: time slope of apparent backscatter; IRC: integrated reflection coefficient; FSIR: frequency slope of integrated 
reflection. 
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Table 3. Pearson's correlation coefficient (r) of parameters and P for Experiment 2. 





AIB (dB) 


FSAB (dB/MHz) 


TSAB (dB/us) 


IRC (dB) 


FSIR (dB/MHz) 


AIB (dB) 
FSAB (dB/MHz) 
TSAB (dB/|!S) 
IRC (dB) 


r = 0.69 e 
r = 0.86 9 
r = -0.56 d 


r = 0.69 e 

r = 0.77 e 
r = 0.72 e 


r = 0.86 9 
r = 0.77 e 

r = -0.50 c 


r = -0.56 d 
r = -0.72 e 
r = -0.50 c 


r = -0.33 a 
r = -0.66 e 
r = -0.37 b 
r = 0.81 f 



a P = 0.24; b P = 0.19; C P = 0.07; d P<0.05; e P<0.01; f P<0.001; 9 P<0.0001. AIB: apparent integrated backscatter; FSAB: frequency 
slope of apparent backscatter; TSAB: time slope of apparent backscatter; IRC: integrated reflection coefficient; FSIR: frequency slope 
of integrated reflection. 



(-4.2 ±0.6 dB/us) to characterize in vitro human trabe- 
cular bone, and reported promising results. In our 
research, the lower values of our AIB parameter, when 
compared with the literature (20), can perhaps be 
explained by the effects of cortical bone, which increases 
the loss of ultrasonic energy by reflection; therefore, less 
energy was transmitted to the inside of the bone, and, 
hence, the energy backscattered from this region was 
decreased. FSAB, on the other hand, had higher (less 
negative) values, when compared with the literature (20), 
indicating that, although the overall energy backscattered 
was smaller, more energy from higher frequencies was 
backscattered. Similar behavior was observed with 
respect to TSAB, where more energy was backscattered 
from deeper structures, thus giving less negative TSAB 
values than those in the literature (20). The low standard 
deviations for parameters AIB, FSAB, and TSAB suggest 
that our method gave statistically similar results at 
different points along the femur; this is clinically 
significant, considering the surface anatomical irregula- 
rities of the bone. 

IRC (13,31,32) represents the amount of reflection 
from a tissue; thus, a denser bone tends to have higher 
values of IRC. Hakulinen et al. (13), using bovine 
trabecular bones in vitro from different regions of the 
femur, found the following values for the parameter IRC: 
-17.9±3.9dB (medial condyle), -20.3±3.8dB (lateral 
condyle), and -27.2±2.0dB (greater trochanter). The 
lower values of IRC observed in our research may be 
because the study had a layer of soft tissue that strongly 



influenced the ultrasonic parameters (33); however, when 
we conducted in vitro experiments (34), we found results 
similar to those in the literature. 

We introduced the parameter FSIR as a new way 
to characterize reflection as a function of frequency (28). 
FSIR showed strong positive correlation to IRC in the 
three experiments, which suggests that this parameter is 
related to the density of cortical bone and is promising 
among the other parameters for monitoring the process 
of bone healing. IRC had a negative correlation to the 
backscatter parameters (except in two cases), indicating 
that it can be an additional parameter for assessing bone 
quality. Furthermore, IRC was negatively correlated to 
parameter AIB. This fact indicates that the method 
proposed in this study succeeded in characterizing the 
surface and inside of the bone. 

The experiments were performed at three different 
periods, with intervals of 10 days, to study the repeat- 
ability of the method. The statistical tests showed no 
difference between parameters AIB, FSAB, TSAB, and 
IRC obtained from the three experiments, but the FSIR 
values obtained in Experiment 1 and Experiment 3 were 
significantly different. This difference might have been 
caused by errors of signal acquisition in a different region 
of interest. Statistical analysis of the correlation between 
the parameters indicated that there was lack of correlation 
in only three cases in Experiment 2; however, overall, the 
experiments suggested that all the parameters have the 
potential to characterize bone. 

The proposed protocol and calculated parameters 



Table 4. Pearson's correlation coefficient (r) of parameters and P for Experiment 3. 





AIB (dB) 


FSAB (dB/MHz) 


TSAB (dB/us) 


IRC (dB) 


FSIR (dB/MHz) 


AIB (dB) 
FSAB (dB/MHz) 
TSAB (dB/us) 
IRC (dB) 


r = 0.96 e 
r = 0.97 e 
r = -0.95 d 


r = 0.96 e 

r = 0.96 e 
r = -0.90 b 


r = 0.97 e 
r = 0.96 e 

r = -0.93 c 


r = -0.95 d 
r = -0.90 b 
r = -0.93 c 


r = -0.86 b 
r = -0.85 a 
r = -0.85 a 
r = 0.90 b 



a P<0.001; b P<0.0001; C P<0. 00001; d P<0. 000001; e P<0. 0000001 . AIB: apparent integrated backscatter; FSAB: frequency slope of 
apparent backscatter; TSAB: time slope of apparent backscatter; IRC: integrated reflection coefficient; FSIR: frequency slope of 
integrated reflection. 
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demonstrate the potential to characterize femur diaphysis 
of rats in vivo by the pulse-echo ultrasonic method 
associated with wave reflection and backscattering. The 
values of the five parameters in rats, as well as the simple 
protocol for signal acquisition, provide supplementary 
predictive data of living human bones. Furthermore, this 
research contributes to the use of rats in vivo in future 
studies of bone characterization to provide data for more 
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